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BACKGROUND  AND  SCOPE 


In  1958  the  General  Engineering  Laboratory  made  a  study  of  propeller 
shaft  thrust  bearing  operation  and  reported  its  findings  in  Reference  1. 
Following  this  study  a  comprehensive  analytical  and  experimental  program 
was  undertaken,  for  the  purpose  of  extending  present  understanding  of  these 
bearings  and  in  order  to  provide  a  body  of  design  information  for  use  in 
bearing  design  and  selection.  This  program,  like  the  preceding  introductory 
study,  was  performed  under  a  contract  awarded  by  the  Bureau  of  Ships  to 
General  Electric  Company's  Medium  Steam  Turbine,  Generator  and  Gear 
Department. 


The  program  was  divided  into  three  phases  as  follows: 


ptt_*ae  I;  Investigate  analytically  the  performance  of  propeller  shaft  thrust 
bearings  using  the  existing  Reynolds -Energy  Method  of  solution 
to  provide  data  necessary  in  design  and  selection  of  these  bearings. 

This  phase  was  completed  and  a  final  report  on  it  was  issued  in 
May  1959- 

Phase  II:  Extend  existing  analytical  techniques  for  propeller  shaft  thrust 
bearings  by  including  a  numerical  method  of  solution  of  the 
Elasticity  Equation.  Review,  and  where  necessary,  modify  the 
design  data  obtained  in  Phase  I,  so  as  to  include  the  effects  of 
pad  distortion  caused  by  pressure  distribution  and  thermal 
gradients. 

This  phase  has  also  been  completed  and  the  following  is  our 
final  report  on  it. 


Phase  III:  Instrument  thrust  bearings  on  two  U.S.  Navy  ships  and  obtain 
experimental  data  on  the  performance  of  these  bearings.  This 
data  is  to  be  obtained  in  tests  carried  out  at  the  time  of  scheduled 
sea  trials.  The  thrust  bearing  performance  measurements 
obtained  in  these  sea  trials  is  to  be  used  for  correlation  with  the 
design  data  obtained  analytically  in  Phases  I  and  II. 

Sea  trials  were  made  during  1959  on  the  U.  S.  S.  Barry,  and  a 
final  report  on  these  tests  is  presently  being  prepared.  The 
second  ship  for  the  sea  trial  tests  has  not  yet  been  selected. 

A  fourth  phase  which  included  the  building  of  a  thrust  bearing  test  stand  was 
contemplated  but  was  not  included  in  the  present  program,  since  the  findings  of 
the  program  could  be  used  to  determine  the  features  of  the  stand. 
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I.  INTRODUCTION: 

The  need  for  greater  accuracy  in  predicting  the  performance  of  tilting 
pad  thrust  bearings  was  discussed  in  References  1  and  2.  In  Reference  2, 
various  performance  calculation  methods  were  described  and  it  was  shown 
that  the  load  capacity  of  thrust  bearing  pads  is  significantly  influenced  by 
their  elastic  deflections,  particularly  in  the  case  of  centrally  pivoted  pads. 

An  approximate  method  of  calculating  the  pad  deflections  and  including  them 
in  the  film  shape  equation  was  worked  at  in  Reference  2  and  it  formed  the  basis 
of  the  design  charts  and  other  calculation  results  given  there.  While  this 
method  of  calculation  gave  results  that  corresponded  more  closely  with  field 
evidence  than  did  results  obtained  by  previous  methods,  it  suffered  from 
several  limitations  the  principal  ones  being: 

1.  Only  deflections  due  to  load  were  considered.  No  allowance  was 
made  for  deflections  due  to  thermal  gradients. 

2.  The  effects  of  cut-outs  and  ring  supports  were  not  considered.  An 
average  pad  thickness  was  used  in  each  calculation  and  point  support 
was  assumed. 

3.  It  was  assumed  that  initially  flat,  pivoted  pads  assume  a  spherically 
crowned  shape  under  load,  the  magnitude  of  the  crown  being 
proportional  to  the  load. 

In  the  study  which  is  reported  here,  a  more  complete  analysis  is  made  in 
which  the  elastic  deflections  of  the  pads  are  calculated  from  numerical 
solution  of  the  Elasticity  Equation.  The  above  three  limitations  of  Reference  2 
are  eliminated  here.  Deflections  due  both  to  load  and  to  thermal  gradients  are 
calculated  and,  in  each  case,  the  actual  pad  geometry,  including  cut-outs  and 
support  ring  is  used  in  the  computations.  Finally,  simple  crowning  of  the  pads 
is  not  assumed;  instead,  the  more  complex  deflected  shape  of  the  pads  is 
described  by  series  expansions  in  r  and  9,  whose  coefficients  are  calculated 
in  each  case  to  satisfy  the  minimum  energy  principle. 

The  analysis  and  calculations  described  in  this  report  had  the  following 
objectives: 

1.  To  derive  a  general  and  accurate  method  for  calculating  the  performance 
of  tilting  pad  bearings  of  any  geometry. 

2.  To  calculate  and  plot  design  charts  for  bearings  with  standard  geometry 
pads  within  the  range  of  sizes  and  speeds  used  in  today's  marine  main 
propulsion  units. 
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3.  To  calculate  the  effect  of  varying  (from  the  standard)  several  of  the 
geometrical  parameters  of  tilting  pads.  This  serves,  at  least,  to 
point  out  the  direction  in  which  these  parameters  should  be  changed 
in  order  to  achieve  more  nearly  optimum  configurations. 

To  meet  these  objectives  an  IBM  704  digital  computer  program  was 
written  which  allows  simultaneous  solution  of  the  Reynolds,  Energy  and 
Elasticity  Equations  for  sector  shaped  pads  of  arbitrary  geometry.  Since  we  are 
concerned  with  tilting  pads,  the  program  includes  a  section  which  calculates 
at  each  operating  condition  the  circumferential  and  radial  pad  inclinations 
which  satisfy  moment  equilibrium. 

This  program  was  then  used  to  calculate  several  bearing  configurations 
at  a  number  of  operating  conditions,  as  follows: 

1.  Four  bearings  with  standard  geometry  pads,  scanning  the  size  range 
25"  O.  D.  to  50"  O.  D.  were  analyzed  at  speeds  up  to  320  RPM.  The 
results  of  these  calculations  were  used  to  prepare  the  design  charts. 

2.  The  ahead  and  astern  bearings  of  the  DD  933  were  analyzed  to  permit 
direct  comparison  of  theoretical  prediction  with  the  field  data  obtained 
under  Phase  III.  (Performance  data  obtained  with  instrumented  thrust 
shoes,  during  the  sea  trials  of  the  DD  933.) 

3.  As  in  Reference  2,  the  31"  O.  D.  bearing  was  selected  to  study  the 
effect  of  varying  bearing  geometry  and  oil  inlet  temperatures. 

Calculations  were  made  varying  the  following  parameters  one  at 

a  time: 


a. )  Pad  thickness 

b. )  Groove  temperature 

c. )  Diameter  ratio  (O.  D. /I.  D. ) 

d.  )  Number  of  pads 

4.  In  order  to  checkout  the  effect  of  the  assumptions  made  in  Reference  2, 
some  additional  calculations  were  made  to  evaluate  the  effects  of 
thermal  distortion  of  the  pads  and  of  the  ring  support  and  of  the  cut¬ 
outs  in  the  pads. 

In  the  succeeding  sections  of  this  report  the  method  of  analysis  and  the 
calculation  results  are  presented  and  discussed. 
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II.  ANALYSIS  AND  METHOD  OF  CALCULATION 


The  following  Assumptions  are  made  in  the  analysis: 

1.  Steady  state  conditions  prevail  in  the  film. 

2.  Lubricant  is  Newtonian  and  incompressible. 

3.  Flow  in  clearance  spaces  is  laminar  and  adiabatic. 

4.  Lubricant  inertia  forces  are  neglected. 

5.  Pressure  and  shear  effects  on  viscosity  are  neglected. 

6.  Variations  of  the  specific  heat  of  the  lubricant  with  pressure 
and  temperature  are  neglected. 

When  the  first  four  assumptions  listed  above  are  valid,  the  hydrodynamic 
pressures  (which  are  generated  in  the  fluid  film  that  separates  the  bearing 
pads  from  the  runner)  may  be  calculated  by  integrating  Equation  1  (Reynolds 
Equation) : 


In  order  to  integrate  Eq.  (1)  for  any  given  geometry  and  speed,  the 
viscosity  and  film  thickness  need  to  be  known  as  functions  of  position.  For 
petroleum  oil,  in  the  range  of  present  day  thrust  bearing  practice,  shear  and 
pressure  effects  on  viscosity  are  very  small  (Ref.  3).  The  local  viscosity  values  are 
then  functions  only  of  the  local  film  temperatures  and  the  viscosity-temperature 
characteristics  of  the  oil.  In  the  present  analysis,  adiabatic  conditions  are 
assumed  and  the  temperature  distribution  is  obtained  by  integrating  Eq.  2 
(Energy  Equation,  Ref.  4)  simultaneously  with  Eq.  1. 


(wr)2  + 


sir[(!fe)  +  (&)  ]  -cpp8*  (nr-njr  ife]  = 

(2) 


The  film  shape,  for  pivoted  pads,  is  determined  by  the  tilt  of  the  pads, 
radially  and  circumferentially  and  by  the  deflections  which  result  from  the 
hydrodynamic  pressures  over  the  pads  and  the  temperature  gradients  across 
the  pads.  In  the  present  analysis,  the  pad  deflections  are  calculated  from  the 
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integration  of  Eq.  3  (Ela.tioity  Equation,  Ref,.  5  ,„d  6). 

.illA 

Et  1 


where: 


Et 

i.  d2  ' 

lit2  r  57  .* 


(3) 


,  E’:  “  ln“gr“’‘‘  '°r  b°,h  ““  *"1  •“«»«  pm...  Local  va lue a 

*"  “"d  “  ,h“  ,h'  '«*'*■  -  <**”««•  m  geometry  are  coh.fd.red. 
The  boundary  conditions  used  with  the  above  equations  are: 

a.  )  Pressures  fall  to  zero  around  the  pad  perimeter. 

b. )  Pad  inlet  oil  is  at  groove  mixing  temperature. 

c- )  Radial  temperature  gradients  are  zero  along  the  inner  and  outer 

circumferences  of  the  pads,  because  of  the  cooling  effect  of  the 
surrounding  oil. 

0. 1  Shea,  ,nd  bendiug  mom.„„  . . .  p,d 

S.no.  pad  deformation,  under  high  i.ad.  „,y  p,oduc.  diverging  flta 

‘  . .  h'"“  br““~»  “»  fll«l  near  ft.  ,d,e.  of 

pad.,  an  additional  condition  ia  impoaed.  Ih, a  .UR,,  fh„  *.  hydrodynamic 
pre.au,..  over  the  pad.  nowhere  .m.„.,  m.n  efmo.pherfc. 

mem  J"'  ‘b°V‘  ‘hr"  *"  .imultaneouaiy,  num.rieaf 

od.  and  an  „,r.„ve  procedure.  Einit,  difference.  „..d  f„,  ,h. 

:dT*t  0',hd  “*'"“‘d’  *-»  Equation.,  ..  ^  „  R.„.  , 

principle,  as  described  in  Appendix  1. 

The  iterative  procedure  used  for  the  simultaneous  solution  of  the 
equations  i*  briefly,  as  follows: 

1-  For  a  given  geometry,  speed  and  minimum  film  thickness,  integrate 

V  ld*  and  Energy  Equations,  using  a  preliminary  estimate  of 
e  elections.  This  integration  is  carried  out  (as  described  in 
efs.  7  and  8  )  and  the  pad  inclinations  are  obtained  as  those  values 
which  yield  a  pressure  profile  whose  center  of  pressure  passes 
through  the  axis  of  the  pivot. 

2-  The  Elasticity  Equation  is  then  integrated  for  the  pressure  and 
temperature  profiles  obtained  in  step  1. 


3.  Reynolds  and  Energy  Equations  are  integrated  once  again,  but 
this  time  for  a  pad  which  has  the  elastic  deflections  calculated  in 
step  2. 

The  above  procedure  is  repeated  until  the  differences  between  two 
successive  iterations  fall  within  a  specified  small  limit. 

From  the  final  values  of  the  film  thickness,  pressure  and  temperature 
profiles  obtained  as  described  above,  the  remaining  items  of  interest  such  as 
total  load  carried  on  the  pads,  oil  flow  and  horsepower  loss  are  calculated  in 
the  same  manner  as  in  Ref.  2,  The  input  data  necessary  for  the  calculations 

described  in  Appendix  2.  The  estimated  accuracy  of  the  calculations  is  discussed 
in  Appendix  3 . 
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III.  RESULTS 

The  method  of  Analysis  outlined  in  the  previous  section  was  used  to 
determine  thrust  bearing  performance  in  the  range  of  sizes,  speeds  and 
configuration  outlined  in  the  introduction  to  (his  report.  The  cases  cal¬ 
culated  are  listed  in  Table  I  and  the  results  are  summarized  in  Tables  II 
through  XVII. 
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IV-  Design  Charts 

The  results  of  the  calculations  made  for  the  bearing*  with  ,  , 

geometry  pad,  (Tah.e,  „  V>  were  p,o„ed  ,  ...  ,  . 

for  ready  use.  et  of  de81gn  charts 

The  average  film  temperature  was  found  from  the  ,  , 

*»  >»  .  Of  sis,  „d  speed.  „  1 J 

plotted  versus  unit  load  in  Figure  1. 

The  maximum  film  temperature  is  olotf.rf  , 

through  10,  for  the  investigated  range  of  beari  Ver8US  8PCed  ^  FigUre<>  2 

made  for  6,  8  and  10  oad  h  ^  81Ze“'  SeParate  charts  were 

o,  and  10  pad  bearings  and  at  three  values  of  • 

thickness  (0.  0014  in. ,  0.  0010  in.  and  0.  0005  i„.  ,  ^  ^ 

-  v",u- — - — -  — > 

for  6.  8,  and  10  pad  bearings  and  for  Jsame  TrlT^slT  ^  ^ 
thickness  (0.  0014  in. ,  0. 0010  in. ,  0.  0005  in. ).  ™nimum  film 

The  hydrodynamic  oil  flow  and  tv...  u 

plotted  in  dimensionless  form  in  F‘  °r8epower  lo»»  ln  the  fluid  film  are 

The  use  of  the  charts  is  illustrated  in  the  following  examples: 

— - *— 

21  1/2  in.  I,n  a  .  lR1  horseP°w«  loss  in  a  43  in.  O.  D.  x 

•  *  pad  thrust  bearing  running  at  280  rpm  ^ 
of  610  psi.  g  280  RPM  and  carrying  a  unit  load 


?avg  at  hmin  “  °-0014  in.  (per  Figure  14) 
Pavg  *  ^min  "  0-001°  in.  (per  Figure  15) 


507  psi 


665  psi 
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Pavg  at  hmin  =  °' 0005  in-  (per  Figure  16) 
Pavg  at  **mjn  =  610  psi  (by  interpolation) 


T  at  h  .  =  0.0014  in.  (per  Figure  5) 

max  min  6  ' 


T  at  h  ,  5  0.0010  in.  (per  Figure  6) 
max  min  r  6  ' 


T  at  h  =  0.0005  in.  (per  Figure  7) 
max  min  6  ' 


^max  at  ^min'*  0*  00114  in.  (by  interpolation) 


^avg  a*  ^10  P8*  ^Per  figure  1) 

^avg  at  Tavg  =  175°F  (per  FiSure  A-20) 

U  =TT(D-L)N  = 
avg 

B  =  1/2  (D-L)  0_  __  = 
avg  TOT 

H  U 

avg  avg 

P  B 

avg  avg 


B  LU 
avg  avg 


(per  Figure  21) 


Q  per  pad 


Total 


6600  HP  (per  Figure  24) 

B  Lp  U 
avg  "avg  avg 

H.  P.  per  pad 

H.  P.  Total 


848  psi 
0.00114  in. 

199°F 

223°F 

286°F 

214°F 

175°F 

2.  2  x  10  Sb  sec/in2 
473  in/ sec 
10.77  in. 

0.  158  x  10'6 

30.  0  x  10‘6 

1.  64  G.P.M. 

13. 1  G.P.M. 

1.0  x  ^ 


5.06  H.  P. 
40.  5  H.  P. 
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V.  Discussion 

The  two  limiting  criteria  of  thrust  bearing  performance  are  the 
maximum  film  temperature  and  the  minimum  film  thickness.  Generally, 
however,  the  main  emphasis  in  the  literature,  as  well  as  in  design  practice, 
has  been  on  film  thickness  with  little  attention  devoted  to  maximum  temperature. 

This  is  because  isothermal  bearing  analysis  was  used  almost  exclusively  in 
the  past.  The  results  of  the  more  complete  study  described  in  this  report, 
in  which  the  effects  of  temperature  (  and  hence  viscosity)  variations  in  the  film 
and  of  the  elastic  deflections  of  the  pads  induced  by  the  thermal  gradients 
were  both  included,  has  convinced  the  authors  that  maximum  temperature  is 
the. governing  limitation  on  bearing  performance.  ThuB,  the  design  charts  of  the 
previous  section  were  prepared  so  that  both  maximum  film  temperature  and 
minimum  film  thickness  can  be  readily  obtained  for  any  operating  point  within 
the  range  studied.  In  the  subsequent  discussion  we  will  attempt  to  show  the 
importance  of  temperature  in  determining  thrust  bearing  performance.  It  will 
also  be  shown  that  when  the  film  and  groove  temperatures  are  known,  either 
analytically  or  experimentally,  the  minimum  film  thickness  also  is  automatically 
defined. 

1.  Maximum  and  Average  Temperatures 

To  illustrate  the  large  temperature  variation  in  the  film,  consider  a 
31  in.  O.D.  thrust  bearing  whose  pad  dimensions  are  given  in  Figure  A-2. 

The  performance  data  for  this  bearing  are  given  in  Table  III  and  the  temperature 
data  has  been  summarized  on  Figure  26.  This  figure  shows  that  while  the 
average  film  temperature  and  the  average  pad  outlet  temperatures  are  comparable, 
they  are  considerably  below  the  maximum  film  temperature,  particularly  at 
high  loads.  In  addition  note  that  the  slope  of  the  maximum  film  temperature  increases 
with  load.  It  is  clear  therefore  that: 

(a)  The  bearing  oil  drain  temperature,  which  is  often  used  to  monitor 
the  condition  of  a  bearing,  is  in  fact  a  very  poor  indicator  of  im¬ 
pending  bearing  failure. 
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b. )  With  a  heavily  loaded  bearing,  even  a  moderate  increase  in  load 
may  raise  the  maximum  temperature  to  a  dangerous  level. 

2-  Influence  of  Groove  Temperature 

A  key  point  of  this  report  is  brought  out  in  Figure  27,  which  shows 
the  relation  between  the  minimum  film  thickness  and  the  maximum  temperature, 
for  a  range  of  groove  temperatures  and  unit  loads.  This  figure  shows  that 
for  the  given  bearing  geometry  and  speed,  the  maximum  temperature  can  be 
used  to  predict  the  minimum  film  thickness  with  an  accuracy  of  better  than 
+  0.0001  in.  If,  in  addition,  the  groove  temperature  is  also  known,  then  the 
minimum  film  thickness  can  be  predicted  even  more  closely.  Thus,  maximum 
temperature  can  be  employed  to  measure  minimum  film  thickness.  This  is 
quite  significant  for  such  temperature  measurements  are,  in  general,  much 
easier  to  make  than  direct  measurements  of  film  thickness  with  presently 
available  indicators.  Also,  based  on  the  above,  their  potential  accuracy  in 
determining  minimum  film  thickness  appear  to  be  superior  to  presently 
available,  direct  indicators  of  film  thickness. 

Fig.  27  also  shows  the  very  important  influence  of  groove  temperature 
on  bearing  performance.  To  illustrate  this  point  note  the  gain  that  can  be 
achieved  by  a  30  F  reduction  in  groove  temperature  when  the  bearing  is 
operating  under  700  psi  load  (data  from  Fig.  27): 


P  -psi 

avg 

T  -  °F 

GR 

h  -  in 

min 

T  -°F 

max 

700 

160 

0.00063 

240 

700 

130 

0.00093 

197 

3.  Influence  of  Thermal  Distortion 

The  influence  of  thermal  distortion  of  the  pads  on  bearing  performance  is 
shown  in  Figs.  28  and  29.  In  these  figures,  the  minimum  film  thickness  and 
maximum  temperature  are  plotted  versus  unit  load  from  the  results  obtained 
using  three  methods  of  analysis: 
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a.  )  Present  analysis  -  complete  solution,  including  the  thermal 

distortion. 

b. )  Present  analysis,  but  neglecting  thermal  distortion. 

c.  )  Simplified  analysis  of  Ref.  2  (In  this  method  thermal  deflections 

were  neglected  and  it  was  assumed  that  the  pads  deform  to  a 
spherical  shape  under  load.  Using  a  simplified  representation 
of  the  pads,  the  radius  of  curvature  of  the  deformed  pad  (R  ), 
was  related  to  the  load  (W),  the  average  pad  thickness  (t)  and  the 
modulus  (E)  by  the  expression: 

R  =2.23  Et3  . ) 

C  W 

The  calculations  were  made  in  all  three  cases  for  the  following  pad 
geometry  and  operating  conditions: 

R  =  15.  5in. 

L  =  6.  75  in. 

t  1  2.2  in.  (uniform  thickness) 

N  =  5.  33  rps 
Oil  =  2190  T 

TCR-  per  plot  of  groove  mixing  temperature  in  Fig.  A-20. 

Note  in  Figures  28  and  29  that  at  least  at  loads  up  to  600  p.i,  the  results 
obtained  using  the  simplified  analysis  of  Ref.  2,  were  reasonably  close  to  those 
obtained  using  the  present  analysis,  but  neglecting  thermal  distortion.  When, 
however,  thermal  distortion  is  introduced  in  the  analysis,  the  situation  is 
drastically  different  at  high  loads.  Thus,  at  unit  loads  greater  than  600  psi 
the  complete  analysis,  including  thermal  deformations,  predicts  much  smaller 
minimum  film  thickness  and  higher  maximum  temperatures.  Thermal  dis¬ 
tortions  are  thus  seen  to  have  a  major  effect  on  bearing  performance  and  must 
be  included  when  designing  bearings  to  operate  under  high  unit  loads. 
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Figure  30  shows  the  film  thickness  contour  lines,  at  a  minimum  film 
thickness  of  0.0004  in.  where  thermal  distortions  are  neglected.  Figure  31 
shows  the  contour  lines  for  the  same  minimum  film  thickness,  but  with  the 
thermal  distortions  included.  Note  from  these  figures  that: 

a. )  The  pad  deformations  are  several  times  the  minimum  film  thickness 

b. )  The  thermal  deflections  are  of  the  same  order  as  those  dirtctly  due 

to  the  hydrodynamic  pressures. 

These  figures  underline  again  the  importance  of  elasticity  in  pivoted  pad 
thrust  bearing  analysis  and  the  need  to  include  the  effects  of  thermal  distortions 
in  such  analysis. 

4.  Influence  of  Geometric  Parameters 

Figures  32  and  33  show,  respectively,  the  influence  of  pad  thickness  and 
ratio  on  bearing  performance  (In  Figure  32,  the  designation  "standard 
pad"  refers  to  the  pad  shown  in  Figure  A-2). 

Note  from  Figure  32  that  at  600  psi  (which  is  on  the  high  side  of  the 
range  of  unit  loadings  for  thrust  bearings  in  marine  use)  the  curves  are 
relatively  flat,  showing  that  thickness  changes  will  not  much  effect  bearing 
performance.  When,  however,  operation  under  higher  loads  is  required, 
pad  thickness  is  seen  to  have  a  greater  effect.  Thus,  for  a  bearing  which  is 
required  to  carry  800  psi,  Figure  32  shows  that  a  30#  increase  in  pad  thick¬ 
ness  would  achieve  a  substantial  increase  in  minimum  film  thickness  and  an 
even  more  important  reduction  in  maximum  temperature 

Figure  33  shows  that  there  is  an  optimum  _L_  ratio,  both  from  the  stand¬ 
points  of  maximum  temperature  and  minimum  *hlm  thickness.  While  the 
optimum  L_  ratio  varies  to  some  extent  with  load,  it  is  in  the  neighborhood 
of  0.  5  within  the  range  of  loads  studied. 

\ 

r 
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5.  Comparison  with  Rigid  Pad 

Thus  far  the  discussion  has  been  limited  to  centrally  pivoted  pads  which 
deform  elastically  under  load,  and  because  of  the  thermal  gradients  generated 
by  the  shear  losses.  The  negative  effects  of  large  deformations  which  occur 
under  high  loads  was  pointed  out.  It  must,  however,  be  remembered  that, 
particularly  for  the  case  of  centrally  pivoted  pads,  some  deformation  is 
desirable  for  adequate  load  capacity.  The  reason  for  this  is  that,  in  order  to 
satisfy  equilibrium  of  moments,  the  resultant  of  the  hydrodynamic  pressures  must 
pass  through  the  pivot.  Thus,  for  the  case  of  centrally  pivoted  pads,  the 
pressure  profiles  should  be  nearly  symmetrical  about  the  radial  centerline. 

This  condition  is  achieved  when  the  active  faces  of  the  pads  exhibit  a  slight 
convexity.  Because  of  the  deflections,  pivoted  pads  will  automatically  achieve 
this  convexity  when  run  under  load. 

In  the  case  of  flat  rigid  pads,  it  is  generally  desirable  to  have  an  off- 
center  pivot,  so  that  the  non-symmetrical  pressure  profiles,  associated  with 
flat  inclined  sliders  can  be  achieved.  A  centrally  pivoted,  rigid,  flat  pad,  on 
the  other  hand,  must  rely  on  the  viscosity  variations  over  its  length,  to  allow 
generation  of  hydrodynamic  pressure  which  satisfy  moment  equilibrium. 

The  tabulation  on  the  next  page  compares,  at  one  operating  point,  the  unit 
loads  and  maximum  temperatures  for  an  elastic  centrally  pivoted  pad  with  those 
for  a  rigid  pad  with  optimum  pivot  and  with  a  central  pivot.  Note  that  the  unit 
loads  carried  by  the  elastic  pad  and  by  the  flat,  rigid  pad  with  optimum  pivot  are 
comparable,  whereas  that  carried  by  the  flat,  rigid,  centrally  pivoted  pad  is 
much  smaller. 
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Comparison  of  Rigid  and  Elastic  Pada 


R  =  15.5  in. 

tgr  "  130°F 

L  =  7.75  in. 

h  =  0,  001  in. 

o 

min 

0T=  38.25° 

oil  =  2190  T 

N  =  5.  33  rpa 

Type  of  Pad  Unit  Loading  Max.  Film 

t _ 


(pai) 

Temp.  (°F) 

! 

Rigid  flat  pad  with 
optimum  pivot 

(r  %  =  51;  6  %  =  61) 

P  P 

660 

179 

,  Rigid  flat  pad,  centrally 

pivoted. 

(r J>  =  51;  0  g  -  so) 

1  P 

370 

204 

1 

Elastic,  initially  flat  pad, 
centrally  pivoted 
(thickness  and  pivot  location 
per  Fig.  A-2) 

642 

194 

t 
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VI.  Summary  and  Prinrip,]  Concision, 

*7“  “’ru“  b“rl"!  a“ly*“  w“  *‘v“  *“  - 

..  temperature  variation.  .„d  ela.tte  deflection.  ]oad 

ZT7  ‘ «...  method  oI  analyeia 
»tde  range  of  epe.de.  load,  a„d  bearing  configuration., 
cwiatton  re.uit.,  ,  of  deeign  _  p„p„ed  for 

“  ‘n  the  ■-“*  '“««•  or  preaent  day  marine 

practice.  Ftta  temperature  ...  .Hewn  be  ,  major 

~rr  b°,h  b““-  -  ~~  * — «  -- 

or  .a  ZT" *•* th**  *•  — -  -c,.-. 

1.  Elastic  deflection.  pad.  .ever,.  „me.  tb.  minimum 

‘  a"Ckn"‘-  thermal  deflection.  „ 

•ho.,  due  directly  ,o  dt.  hydrodynamic  p,..„„  ,o,d,„g.  Th.„ 
deflection,  eerion.l,  affect  th,  performance  of  th.  bearing  ,nd 
non.,  be  included  in  the  an.Iy.i.,  Analy.,.  „hieh  no,,,.,,  ^ 
.la.“e,.y  M  the  effect,  o,  thermal  gradient,  give.  .,ro„.o„. 
and  optimise  picture  o,  thru.,  bearing  operation  high 
For  operation  under  high  load.,  improvement  in  p. rform.nc,  can 
•  gamed  by  minimiaing  deflection.  ,„d  by  low,„„g  groove 
mixing  temperature. 

2';i—  r"'h'd  -  “>«  dra.tic.lly  different 

b<"h  ""  *V'r*««  ,ita  r.hrperatur.  and  the  average  drain 
temperamre.  The  elope  of  the  maximum  temperature  ve.  „„i,  lo,d 

. . .  *'  h“h  >“de.  Highly  loaded  bearing,  of 

I"  ^  h‘V"  “”te  —  -  —T  ror  overload 
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3.  Accurate  determination  of  the  minimum  film  thickness  in  a 

bearing  can  theoretically  be  made  by  temperature  measurements 
alone.  It  is  desirable  for  this  purpose  to  measure  both  the  maximum 
film  temperature  and  the  groove  temperature.  This  method  offers 
greater  accuracy  than  present  direct  film  thickness  indicators. 
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VII.  Recommendationa; 

The  authors  feel  that  future  thrust  bearing  investigations  should  be 
aimed  at: 

1.  Modification  of  the  Energy  Equation,  to  include  heat  conduction. 

2.  Analysis  of  thrust  bearings  under  transient  conditions  in  the  fluid  film. 

3.  Analysis  of  thrust  bearings  under  turbulent  conditions  in  the  fluid  film. 

4.  Meticulous  and  well  instrumented  experimental  study  of  tilting  pad 
thrust  bearing  performance.  It  is  realized  that  such  experimental 
work  would  undoubtedly  be  both  difficult  and  expensive.  It  is,  however, 
essential  in  order  to  test  the  validity  of  analytical  work  and  present 
understanding  of  bearing  performance. 

5.  Utilize  the  analytical  and  experimental  results  to  optimize  the  dimensions 
of  bearing  designs  in  present  use  and  to  generate  new  designs  capable 

of  much  improved  load  carrying  capacity. 
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APPENDIX  I 


Method  of  Solution  of  Elasticity  Equation 

1.  Assumed  Conditions 

Aside  from  the  usual  assumptions  of  elasticity  theory,  the  following  have 
been  used  in  the  numerical  solution  of  Eq..  3r.  . 

a.  Only  deflection  due  to  bending  is  considered;  the  contributions  due  to 
shear  deformation  and  to  membrane  and  crushing  stresses  are  omitted. 

b.  The  support  plate,  or  button,  is  circular,  of  uniform  thickness,  and  at 
uniform  temperature. 

c.  The  support  point  is  at  the  center  of  the  support  plate. 

d.  Load  is  transmitted  between  the  support  plate  and  the  sector  pad 
through  line  contact  around  the  outer  edge  of  the  plate;  the  two  parts 
are  in  contact  around  the  entire  circumference. 

e.  No  moment  is  transmitted  through  the  plate-to-pad  joint. 

f.  Sharp  discontinuities  in  the  pad  thickness  are  averaged  out  over 
finite  area  subdivisions;  and  local  stress  concentrations  are  ignored. 

g.  Temperature  of  the  bottom  of  the  pad  is  assumed  to  vary  linearly 
with  radial  and  circumferential  locations. 

h.  Temperature  of  the  top  of  the  pad  is  assumed  the  same  as  the  lubricant 

film  temperature. 

i.  No  attempt  is  made  to  consider  actual  heat  flow  distribution;  the 
temperature,  in  effect,  is  assumed  to  act  as  if  it  varies  linearly 
through  the  pad  thickness  at  every  point. 

2.  Coordinate  Systems 

The  three  sets  of  coordinates  used  are  shown  in  Fig.  A-l.  Axis-of- 

rotation  coordinates  (  r  ,  6  )  ar«  used  £or  "ubdividing  the  pad  and  for  calculations 

c 

based  on  these  subdivisions,  including  all  integrations  over  the  pad  surface  and 
boundary.  Point-of-support  coordinates  (r^  M  are  used  for  description  of  the 
deformation  modes  of  both  the  pad  and  the  support  plate.  Rectangular 
coordinates  (X,  Y)  are  used  for  load  moment  and  center  of  pressure  eccentricity 
calculations.  The  location  of  the  support  point  (  rp,  6p).  the  circle  of  contact 
between  pad  and  support  plate  (rc  =  tq),  and  the  boundary  coordinates  of  the 
pad  are  also  shown. 


V 
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3.  Energy  of  Bending 

The  selection  of  deformation  mode  functions  is  critical;  it  was  made 
in  this  case  with  consideration  of  the  symmetry  of  the  support  pad  and  of  the 
unimportance  of  the  absolute  level  of  the  system.  Consequently,  the  reference 
for  deflection  calculation  is  taken  as  the  average  location  of  the  edge  of  the 
support  plate;  and  the  shape  of  the  pad  is  written 


where  A.  is  a  coefficient  of  length  dimension,  and  f.  is  a  dimensionless  mode 
shape . 


The  actual  deflection  mode  shapes  for  the  pad  were  chosen  from  among 
the  axially  symmetric  and  wave  deflection  modes  arising  in  the  solution  in 
polar  coordinates  of  the  plate  bending  equation 


for  uniform  plate  thickness,  and  axially  symmetric  or  wave  distributions  of 
load  q.  These  modes  were  classified  as 

Symmetric : 


•  =  0,  2,  3,  4,  . . . 


(A-3) 


Non-symmetric  Type  I: 
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Solution*  involving  In  *  are  not  included  i„  ,hi.  .l,ho„,h  ,hey 
b.  u..d  i,  dealred.  The  Type  I  „o„-.ymm,lric  tt,m,  „  be  ^ 

one.  involving  non-eero  d.il.ction  th,  .upper.  pl„e  a„d  „ 

give  the  only  coupling  with  that  plate . 

The  boundary  condition  on  the  support  plate  is,  therefore,  of  the  form 

'sin  n 

(A-6) 


:os  n  0 


eorre.ponding  ,o  any  Type  1  term..  ror  .  uniform  ^  ^ 

load  giving  .„ch  an  edge  deflection,  bending  nation  „„  be  aolved  with 
this  result: 


V’  A  n(n-l)(l-v)  /  r  \  n  I  ,  ?| 

*  4_  JilrnTv)  \f)  m 

'  '  I  \  o  /(  ^cos  n  0 


(A-7) 


where  v  is  Poisson's  ratio. 


in  polar  coordinate,.  ,he  bending  energy  of  a  pl„e  „s„  ^ 

written  as  the  following  surface  integral: 


v=  /  y  J(g  +  e)2  -2  (1-v)  ge  +  2  (1-v)  k2 


dS 


where 


Et 


12  (1  -  v) 

i.  the  stiffness  in  term,  of  Young's  Modulus  E  and  plate  thickness  t;  and 


(A-8) 


(A-9) 


g  = 


1  dw 

7  37 

*  h 

r 

d2W 

ae2 

(A-10) 

S2W 

dr2 

(A-ll) 

d  / 1 

dw  1 

§7(7 

3?  J 

(A-12) 
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are  the  deflection  derivatives.  In  the  case  of  the  support  plate,  the 
deflection  (A-7)  can  be  substituted  and  gives  the  bending  energy 

v  -  v  a2  n2(n2-l)(l-v){3+ V) 

SP  i  1  2r  2  (2n+l  +  v) - 

o 

The  variation,  which  is  of  interest  in  the  Ritz  method,  is 


$V. 


SP 


Z  A  6A  —  n2(n2-l)(l-v)(3tv) 

j  i  i  2  (2n  + 1  +  v) 


(A-13) 


(A-W) 


In  the  case  of  the  bearing  pad,  the  nature  of  the  cutouts  and  the  shape 
of  the  pad  make  direct  integration  for  the  energy  impractical.  The  result  is, 
therefore,  approximated  by  computing  an  average  or  mean  value  of  the  integrand 
for  each  of  the  subdivisions  of  surface  used  in  the  lubrication  calculation;  and 
these  are  summed  using  the  weight  factors  which  were  developed  for  the 
integration  of  pressure  to  give  load.  These  weight  factors  are  not  quite 
appropriate  for  the  bending  energy,  but  seem  within  the  accuracy  of  other 
approximations.  Actually,  it  is  convenient  to  program  the  calculation  of  the 
variation  of  the  bending  energy  directly;  symbolically, 


6V'.«  Z 
i.J 


!ij  Ai  SAj 


(A-15) 


The  coefficients  I  correspond  to  integrals  of  the  form 


‘IJ 


<gi  +  "j)  (gj  +  -  (1  -  v){gi  +  gj  et)  +  2  (1  -  v)  kj  k^  j  dS  (A-16) 


In  this  case,  both  subscript  i  and  subscript  j  refer  to  the  individual  deflection 
modes  and  can  range  over  all  of  the  modes  which  are  under  consideration  for 
the  pad.  In  general,  the  cross-coupling  terms  will  not  cancel,  as  they  do  in 
the  case  of  the  support  plate;  so  that  the  terms  ofT^  can  form  a  complete 
matrix  of  order  equal  to  the  number  of  assumed  modes.  The  Ty  matrix  is 
obviously  symmetric,  however,  limiting  the  number  of  combinations  that  have 
to  be  calculated  to  that  extent. 
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The  variation  of  the  total  bending  energy  is  the  sum  of  the  pad  and 
•  upport  plate  values: 


6V  =  6Vp 


+  6V_ 


SP 


(A-17) 


The  result  of  this  summation  is  a  matrix  1  of  coefficients  for  the  parameters 
A.  6 Ay  which  differs  from  the  pad  matrix  1  only  in  those  main  diagonal  terms 
corresponding  to  type  I  wave  modes  coupling  with  the  support  plate.  The 
resulting  matrix  I.j  is,  therefore,  also  symmetric.  In  the  present  case,  this 
represents  the  total  elastic  energy  variation;  and  since  it  depends  only  upon 
the  assumed  deflection  modes  and  the  pad  and  plate  geometry,  it  can  be 
computed  initially  before  starting  any  of  the  iterative  stages  of  the  bearing 
program.  Moreover,  for  the  point-supported  bearing  it  turns  out  that  I.  is 
the  complete  matrix  of  coefficients  of  the  linear  equations  which  have  to  be 
solved  for  the  deformation  parameters  A.;  consequently,  the  initial  calculations 
can  also  include  the  inversion  of  Iy  to  the  solution  matrix  J  ,  giving  a  further 
reduction  in  the  amount  of  calculation  required  in  the  iteration  loops. 


4.  Load  Energy  Variation 

The  energy  variation  corresponding  to  the  lubricant  pressure  loading  is 
given  by  the  variation  of 


Wp  =  -  /  P  X  Aj  fj  <rc.  ♦c)  dS  (A_18] 

where  P  is  the  pressure  differential  across  the  pad,  and  is  equal  to  the  pressure 
i.omputed  from  the  solution  of  the  lubrication  equations.  The  support  plate  is 
assumed  free  from  pressure  differentials.  Since  the  pressure  is  known  only 
at  discrete  points,  the  centers  of  the  mesh  areas,  numerical  integration  is 
again  indicated.  This  is  carried  out  for  each  deflection  mode  in  the  same  way 
as  the  integration  of  pressure  to  obtain  total  load.  The  result  gives  the 
constants  K ^  of  the  expression 


”  =  f 


VAj 


This  is  the  only  mechanical  loading  that  needs  to  be  considered. 


(A-19) 
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There  is,  however,  a  contribution  to  the  deflection  due  to  the  temperature 
variation  through  the  pad  thickneaa;  any  poaaible  similar  effect  in  the  support 
plate  has  been  neglected.  Rather  than  attempt  a  rigorous  anlaysis  of  this  effect, 
the  assumptions  were  made  with  the  objective  of  conceptually  replacing  the 
temperature  effect  by  an  equivalent  mechanical  loading,  the  energy  variation 
then  being  relatively  simple  to  write.  Thus  the  listed  assumptions  effectively 
describe  the  temperature  field  as  varying  linearly  through  the  pad  thickness  at 
every  point.  Deformation  in  the  plane  of  the  pad  could  be  prevented  under  these 
circumstances  by  developing  an  internal  bending  moment,  expressed  in 
rectangular  coordinates,  of 


M  *  hi  =  D(l+ v) 
x  y  t 

where  «<  is  the  thermal  expansion  coefficient, 
length  in  the  appropriate  directions. 


(A-20) 


These  moments  are  per  unit 


The  development  of  such  internal  moments  requires,  in  effect,  a  distributed 


moment  loading  of 

mx=*(1+v)«  It 


(A-21) 


m  =  -  (1+  v)  «< 

y 


d 

57 


(A-22) 


each  of  these  being  expressed  per  unit  surface  area;  and,  in  addition, 
concentrated  edge  moments  around  the  free  boundary  given  by  expression 
(A-20).  The  total  effect  of  the  temperature  distribution,  and  of  the  fore¬ 
going  distributed  and  concentrated  moment  loads,  is  to  maintain  the  pad 
and  at  the  same  time  in  moment  equilibrium. 


The  deflection  due  to  temperature  can,  therefore,  be  represented  by  the 
effect  of  applying  moment  loading  equal  and  opposite  to  that  just  described. 
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The  potential  energy  associated  with  such  a  load  can  be  written  ae 

WTX/  <-mx  5T-my  dS-  /  Mi  sJ  d&  <A'23> 

the  line  integral  being  taken  around  the  boundary  of  the  pad.  Again,  as  a 
practical  matter,  the  integrals  are  approximated  by  appropriate  summations, 
and  the  energy  variation  is  calculated  directly.  Allowance  for  the  description 
of  temperature  variation  on  the  bottom  of  the  pad  is  made  by  writing  the 
temperature  differential  as 

iT  ■  T  -{TGK+  "8  •  9  ♦  «,  -5^  ■  ['■-  <*-«]]  «*-“> 

where  T_^is  the  oil  inlet  edge  temperature;  and  q  ,  , 

allow  for  thefirst  power  variation  along  the  radial  and  circumferential  directions. 

The  (x,  y)  coordinates  indicated  in  (A-20)  through  (A-23)  need  not  be 
rigorously  identified  with  a  fixed  coordinate  system,  as  their  function  is 
ultimately  only  to  give  means  of  computing  energy  density,  which  is  then 
integrated  by  the  summation  method.  In  the  program,  they  are  interpreted 
as  being  the  local  radial  and  tangential  directions  at  each  mesh  center  or 
boundary  point,  since  these  were  the  best  directions  for  carrying  out  the 
necessary  numerical  differentiations. 


The  result  of  the  foregoing  temperature  calculations  is  another  set 
of  load  potential  energy  terms  corresponding  to  6WT;  the  total 


6W  a  6W  +  6Wm 
P  T 


(A-25) 


comprises  a  set  of  coefficients  K  for  the  parameters  <5A  . 

*  J 


5.  Deflection  Equations 

By  the  minimum  or  stationary  energy  principle,  the  equilibrium  state 
is  approximated  when 


6{V  +  W)  =  0 
or 


(A-26) 


2  I  A  6 A  +  2  K  6A  =  0 
i,j  ‘J  1  J  j  j  j 


(A-27) 
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Since  the  variation*  of  the  A'*  are  arbitrary,  thie  i*  equivalent  to  the  *et 
of  linear  equation*  represented  by  the  matrix  equation 


M  H  *  N  =° 

The  solution  i* 

H  •  H  "  H  ' 


.,]  ■  ■  ['./I  i'.. 


(A -28) 


(A-29) 


(A-30) 


The  load  energy  vector  K  i*  calculated  at  each  appropriate  point  in  the  program, 
and  multiplied  by  the  initially  calculated  inverse  matrix  J.j  to  get  the  new  set 
of  deflection  coefficients  A^. 


6.  Functions  Used 

In  the  program  development,  and  in  the  calculations  made  so  far,  six 
deflection  modes  have  been  used.  They  are  as  follows: 


Symmetric: 

fl=  I 


2 


1 


f 


2 


-  1 


Non-symmetric  Type  I: 


Non-symmetric  Type  II: 


sin  2  ♦ 

c 


(A-31) 
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The  choice  was  made  somewhat  arbitrarily,  to  include  all  the  types  of 
functions,  and  to  match  to  as  good  an  approximation  as  possible  what  is 
known  and  suspected  about  pad  deformation.  The  validity  of  the  choice 
has  not  been  tested,  but  operation  of  the  method  appears  to  be  successful 
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APPENDIX  II 

input  data 


1-  Pad  Geometry 

As  described  in  Appendix  I,  the  programed  solution  of  the 
Elasticity  Equation  takes  account  of  the  thickness  variations  in  each  pad 
and  of  the  fact  that  the  pads  are  ring  supported  on  circular  disks.  In  the 
calculations,  the  pad  is  divided  into  a  13  x  13  field  and  its  thickness  profile 
is  represented  by  the  set  of  thicknesses  in  each  mesh.  A  typical  pad  is 
illustrated  in  Figure  A-2  which  shows  the  information  on  pad  geometry  that 
is  introduced  as  input  data  in  each  calculation.  This  data  is  shown  in 

Figures  A-3  through  A-18  for  all  the  geometries  analysed  in  the  work  reported 
here. 

2-  Pad  Inlet  Temperature 

The  oil  temperature  in  the  grooves  between  pads  is,  in  general, 
considerably  higher  than  that  of  the  housing  inlet  pads.  This  is  due  to  the 
mixing  in  each  groove  of  cold  inlet  oil  with  hot  oil  discharged  from  the  pre¬ 
ceding  pad.  The  groove  temperature  is  used  in  the  analysis  as  a  boundary 
condition  of  the  Energy  Equation  (Eq,  2). 

In  the  calculations  described  in  this  report  the  experimentally 
obtained  plot  of  groove  mixing  temperature  vs.  unit  load  that  was  first  used 
in  Reference  2,  was  used  again  to  determine  the  pad  inlet  temperature  in 
each  calculation.  This  curve, which  is  reproduced  here  is  Figure  A-19,  was 
used  in  all  the  cases  calculated,  except  for  casses  110  through  123  (Table  XIIl) 

In  these  cases,  which  were  run  for  the  purpose  of  studying  the 
effect  of  varying  the  groove  temperature,  curves  parallel  to  that  of  Figure  A-19, 
but  displaced  15F  above  and  below  it  were  used. 
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3-  Bottom  Fact  Temperatures  of  the  Pad 

As  shown  by  the  Elasticity  Equation  (Eq.  3).  the  thermal  dis¬ 
tortion  of  a  thrust  pad  is  a  function  of  A  T,  the  axial  temperature  difference 
between  the  oil  film  and  the  bottom  face  of  the  pad  (See  Eq.  A-24).  In  the 
calculations,  the  film  temperatures  were  obtained  from  the  solution  of  the 
Energy  Equation,  while  the  bottom  face  tempe  ratures  were  assumed  to  be 
at  groove  tempe  rature  along  the  inlet  edge  and  to  rise  at  the  rate  of  40  F 
per  radian,  proceeding  towards  the  trailing  edge.  This  gradient  was  ob¬ 
tained  from  the  data  gathered  during  trials  of  the  U.  S.  S.  Barry.  In  these 
trials,  the  back  face  temperatures  of  the  31  in.  O.  D.  bearing  pads  were 
measured  at  several  points  between  the  leading  and  trailing  edges. 

4.  Oil  Properties 

The  lubricant  properties  that  are  used  in  the  calculations  are: 
specific  heat,  mass  density  and  viscosity-temperature  relation.  In  the 
present  calculations,  the  lubricant  properties  used  were  those  of  2190  T 
oil  and  these  are  shown  in  Figure  A-20. 
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APPENjCIX  ITT 

estimate  cf  errors 


A  13  «  ,3  ...  integr„ioii  of  Eqj 

3-  Thi’  W"  *h'  “■“■*  “>*>  ““M  be  „.,d  to  ft,  „ 

”■  *“V  Pr*v‘°",1>'  by  the  .uftor.  or  reported  ft. 

«  ft™,  beta,  C.mp,ri.„„  wa.  n>de  „ 

With  a  calculation  made  using  an  11  *  n  me8h.  Xhe  re  , 

con.  me  resulting  pressure  and 

temperature  fields  agreed  within  1 %. 

Since  an  iterative  method  was  used  in  the  calculations,  limits 

were  set  to  determine  when  the  final  solution  was  achieved.  These  were 
as  follows: 

’  Successive  Heretic,  of  ft,  pre..„„  field  k.d  „iftin 

b.  I  Sucre.., v.  iteration,  of  ele.tic  defies, i„.  h,d  to  ,,p.„  withto 

3  x  10  in. 

c.  )  The  resultant  of  the  oressur.  u-j  * 

P  sure  field  had  to  coincide  with  the 

pivot  within  0.2  x  10 ' 3  in. 

The  authors  feel  that  these  limits  were  sufficient  to  insure  that: 

a. )  Errors  in  calculated  maximum  temperatures  5°F. 

b. )  Errors  in  calculated  minimum  film  thickness  ^0.0001  in. 

much  la  With  regard  t0  tHe  CalCUUted  011  fl°W  3nd  h0rSeP°Wer  *>••■  however, 
Z  J  err°rS  arC  intr0dUCCd  ^  ^  ne—V  numerical  calculation  of 

the  o Tre  gradiCnt8  41  ^  ed8CS  (8Ce  2)’  Ba8Cd  «  ation  of 

the  point  scatter  obtained  from  a  large  number  of  calculations,  it  is  felt  that 

“  ^  Ca8e8’  the  CalCUlatBd  0il  -  Action  losses  may  be  in  error’ 

by  as  much  as  10%.  y  error 
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TABLE  I 


Case  No. 


1 

through 

18 

19 

through 

36 

37 

through 

54 

55 

through 

81 

82 

through 

87 


88 

through 

90 

91 

through 

96 

97 

through 

102 

103 

through 

105 


Bearing  Size  Pad  Speed  Minimum  Film  Calculation  Comments 
Size O^.D. jcI..DAngle  R.P.M.  Thickness  Results  in 

(inches)  (Deg.)  (inches)  Table 


25  x  12  1/2 

51 

160 

0. 0014 

II  ~ 

38.  25 

and 

0. 0007 

and 

320 

and 

30.  6 

0. 0004 

31  x  15  1/2 

51 

160 

0.  0014 

III 

38.25 

and 

0. 0007 

and 

320 

and 

Standa  rd 

30.6 

0. 0004 

Geometry  Bearings 

41  x  20  1/2 

51 

160 

0.  0015 

IV 

38.25 

and 

0. 0008 

and 

320 

and 

30.6 

0.0005 

50  x  25 

51 

100 

0.  0015 

V 

38.  25 

170 

0. 0008 

— 

and 

and 

and 

30.6 

320 

0.  0005 

31  x  16  1/2 

38.25 

160 

0.  0014 

VI 

Ahead  Bearing 

and 

0.0007 

of  DD933.  Calcu- 

320 

and 

lations  made  for 

0. 0005 

comparison  with 

sea  trial  data. 

26  x  17  1/2 

29 

160 

0.  0018 

VII 

Astern  Bearing  of 

0.  0010 

DD933.  Calculations 

and 

made  for  comparison 

0.0004 

with  sea  trial  data. 

31  x  16  1/2 

38.25 

320 

0.  0014 

VIII 

0. 0007 

and 

0. 0004 

Calculations  made 

26  x  17  1/2 

29 

160 

0.  0018 

IX 

»  to  study  influence 

0. 0010 

of  thermal  deforma- 

and 

tion  of  the  pads 

0.0004 

50  x  25 

30.  6 

170 

0.  0015 

X 

0. 0010 

and 

0.  0004 

TABLE  I  (Contd. ) 


Case  No. 

Bearing  Size 
Size  0.  D.  xl. 
(inches) 

Pad 

D. Angle 
(Deg. ) 

Speed 

R.P.M. 

Minimum  Film 

Thickness 

(inches) 

Calculation 
Results  in 
Table 

Comments 

106 

31  x  15  1/2 

38.25 

160 

0.0014 

XI  and 

Calculations 

through 

and 

0.0007 

XII 

made  to  study 

117 

320 

and 

influence  of 

0.0004 

pad  thickness 

118 

31  x  15  1/2 

38.  25 

320 

0.0014 

XIII 

Calculations 

through 

0.0007 

made  to  study 

123 

and 

influence  of 

0.0004 

groove  temp¬ 

erature 

124 

31  x  21 

38.  25 

320 

0.0014 

XIV 

through 

0.0007 

126 

and 

0.0004 

127 

31  x  10 
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0.  0014 

XV 

through 

0.0007 

Calculations 

129 

and 

made  to  study 

0.0004 

influence  of 

130 

31  x  21 

38.25 
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0.0014 

XVI 

pad  geometry 

through 

0.0007 

132 

and 

0.0004 

133 

31  x  15  1/2 

38.25 
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0.0014 

XVII 

through 

0.0007 

135, 

and 

0.0004 
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6,  8  and  10  PAD  BEARINGS 
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PROPERTIES  OF  2190T  OIL 


SYMBOLS 


A, A. 

Deflection  coefficients  (see  Eqs.  A-l,  A-18) 

inches 

B 

avg 

Average  circumferential  pad  length  «  1/2  <D’L)eTOT 

inches 

C 

P 

Specific  heat 

BTU/lb.x°F 

D 

Flexural  rigidity  «  Et^ 

12(l-v2) 

lb.  in. 

E 

Elastic  modulus 

psi 

e 

Deflection  derivative  (see  Eq.  A- 11) 

inches’1 

f 

Deflection  functions  (see  Eqs.  A-l,  A-18,  A-31) 

8 

Deflection  derivative  (see  Eq.  A-10) 

inches"1 

h 

Film  thickness 

inches 

hmin 

Minimum  film  thickness 

inche  s 

*1.  J,Ti.J 

Set  of  deflection  coefficients  comprising  matrices 

J 

Mechanical  equivalent  of  thermal  energy  («9339  In.  lba/BTU) 

Inverae  of  1  .matrix 

J 

KJ,lrj 

Set  of  deflection  coefficient,  comprising  matrices 

k 

Deflection  derivative  (tee  Eq.  A-12) 

inches  * 

L 

Radial  length  of  pad 

inches 

M  ,M 
x  y 

Edge  bending  moments  per  unit  length 

in.  lbs/in 

m 

Deflection  modes  (tee  e.g.  Eq.  A-4) 

m  ,  m 
x  y 

Internal  bending  momenta  per  unit  length 

in.  Ibs/in 

me 

circumferential  inclination  of  pad 

inche  e/inch 

m 

r 

radial  Inclination  of  pad 

inche e/inch 

N 

Speed  of  rotation 

RPS 

n 

Deflection  modes  (eee  e.g.  Eq.  A-4) 

P 

Pressure  differential  across  pad 

psi 

P 

Unit  loading 

psi 

p 

avg 

Average  unit  loading 

psi 

Q 

Hydrodynamic  oil  flow  per  pad 

gpm 

qr 

Surface  temperature  gradient,  radially 

°F/inch 

qS 

Surface  temperature  gradient,  circumferentially 

°F/  radian 

R 

Outeide  radius  of  pad 

inches 

SYMBOLS  -  Cont'd. 


o 

o< 

8 

AT 

0 

0 

P 


®TOT 


v 

♦ 

c 

fa) 


Radius  of  support  ring  inches 

Radial  coordinate  inches 

Radial  coordinate  (see  Fig.  A-l)  inches 

Radius  of  support  pad  (see  Fig.  A-])  inches 


Radial  coordinate  of  pivot  inches 

Radial  coordinate  of  pivot  *  100  jr  -  (R-L)]  /L  inches 

Deflection  modes  (see  e.g.  Eq.  A-3) 

Temperature  °F 

Average  film  temperature  °F 

Groove  mixing  temperature  °F 

Maximum  film  temperature  °F 


Pad  thickness  inches 

Thickness  of  support  ring  inches 

Average  surface  velocity  «TT(D-L)  N  inches 


Energy  of  bending 

Energy  of  bending  of  support  plate 

Energy  of  bending  of  pad 

Potential  energy  associated  with  load 
Potential  energy  associated  with  pressure  loading 


lb.  inches 
lb.  inches 
lb.  inches 

lb.  inches 
lb.  inches 


Potential  energy  associated  with  equivalent  tempera¬ 
ture  loading 
Elastic  deflection 
Deflection  of  rim  of  support  plate 

linear  coefficient  of  thermal  expansion 
Increment 

Temperature  difference  across  pad 
Angular  coordinate 
Angular  coordinate  of  pivot 


lb.  inches 

inches 

inches 

i/°r 

°r 

radians 

radians 


Angular  coordinate  of  pivot  ■  100  (0^/ 0T) 
Pad  subtended  angle 

Absolute  viscosity 
Poisson1  s  Ratio 

Angular  coordinate  (see  Fig.  A-l) 


radians 
radians 
lb.  sec/ in* 

radians 


Angular  velocity 


radians/ sec. 
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'  thrue*.  brgs.  is  presented  in  which  the  temperature  dlstri- 
;  but  ion  in  the  fluid  film,  &  the  pressure  &  temperature  in¬ 
duced  pad  deflections  are  considered.  The  method  consists 
of  simultaneous,  numerical  integration  of  the  Reynolds, 
Energy,  6>  Elasticity  Eqs.  Calculation  results  and  charts 
are  presented  covering  the  size  &  speed  ranges  of  present 
propeller  shaft  practice.  The  Influence  of  temperature  and 
the  effect  of  varying  pad  geometry  are  discussed. 
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CONCLUSIONS 

1.  Elastic  deflections  of  tilting  pads  are  several  times 
the  minimum  film  thickness.  Thermal  deflections  are  of 
the  same  order  as  those  due  to  load.  These  deflections 
seriously  affect  bearing  performance  under  heavy  loads  and 
should,  therefore,  be  Included  in  the  analysis. 

2.  The  maximum  temperature  reached  in  the  lubricant  film 
differs  drastically  from  the  average  film  and  the  average 
outlet  temperatures.  The  slope  of  the  maximum  film  temp¬ 
erature  vs.  unit  load  curve,  rises  sharply  at  high  loads 
so  that  heavily  loaded  bearings  have  little  isargln  of 
safety. 

3.  Accurate  determination  of  minimum  film  thickness  is 
theoretically  possible  from  measurements  of  the  maximum 
film  and  groove  temperatures.  This  method  offers  greater 
accuracy  than  present  direct  film  thickness  indicators. 
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